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is discussed fto« three different aspects. First, the 
method is presented; next, the flom chart is traced; and 
finally, the required input is examined. The accuracy of 
the algorithm is evaluated by comparing its predictions of 
reguired cooling flow with those of major engine 
manufacturers. 
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ANALYSIS 

Turbine Cooling Method 

The quantity of required cooling flow and the 
corresponding decrease in stage efficiency are calculated 
ior each row of airfoils throughout the turbine. These 
values are then used to obtain the guantity of conpressor 
bleed flow required to cool the turbine and the decrease in 
cooled-turbine efficiency caused by cooling air injection 
into the gas stream. The calculations depend on both the 
type of cooling configuration and the value of cooling 
ef fectiveness. 


Cooling confiqurat i 
configurations in this 
table I. In this table, 
used to identify the co 
2. Columns 3 and i are 
configuration for each r 
required cooling flow 
that flow requirsd by 
blade. For example, a 
requires only 0.8 of t 
coverage film cooled tu 
cooling effectiveness, 
decrease in stage effic 
required to cool the sta 


,2S* “ The allowable cooling 

cooling algorithm are described in 
column 1 contains NFACT(I) which is 
oling conf iguration shown in column 
aids for choosing the proper cooling 
ow of airfoils. Column 5 lists the 
for each configuration relative to 
a full coverage film cooled turbine 
transpiration cooled turbine blade 
he cooling flow reguired by a full 
rbine blade for a given value of 
Columns 6 and 7 list the percent 
iency for each percent of core flow 
tor(s) and rotor (s) , respectively. 


The values of relative cooling flow, FACTOR, which 
appear in column 5 of table I come from both analytical and 
experimental sources. Reference 3 presents an analysis 
relating the required flow for convection, advanced 
convection, full coverage film, and transpiration cooling. 
These values are sensitive to gas, coolant, and metal 
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teaperdtures and to dirfoil 
purpose of this investiq-itio n, 
cooLiny flow are sufficient. 
FACTOR of 2.0, 1.5, 1.0, and 

to coolinq conf iyurat ions 1, 3 


geometry. However, for the 
itoainal values of relative 
Froa reference 3 values of 
0.8 are obtained correspond i ng 
, 8, and 9, respecti vel y. 


. figure 8 of reference 4, a convection cooled 

turbine blaae with a 0.023 centimeter coating requires only 
three-fourths of the cooling flow of that of an uncoated 
convection cooled blade. So the value of FACTOR of 1.3 for 
a coated convectio.. cooled blade can be obtained from the 
value of 2.0 for a convection cooled blade. 


Values of FACTOR for cooling configurations 4, 5, and 6 
are 1,3, 1.2, and 1.1, respectively. These values have been 
obtained by interpolating between values of FACTOR for 
cooling configurations 3 and 8, The value of FACTOR for 
cooling configuration 7 was obtained by using a weighted 
average of FACTOR values for cooling configurations 3 and 9. 

Values of the dimensionless loss in stage efficiency, 
DBLV and DELM, which appear in columns 6 and 7 of table I, 
were obtained in a personal communication from Thomas p 1 

•loffitt, a member of the Turbine Branch at lewis Research 
Center. 

iZ£S£SS« ~ Cooling effectiveness, PHI 
IS defined as the ratio of the difference between the hot 
gas temperature and the allowable bulk metal temperature to 
the difference between the hot gas temrerature and the 
compressor bleed temperature. 


PHI = (T - TM)/(T - TC) 


The hot gas temperature entering a given ruv of 
airfoils is the average combustor exhaust tempe'*atiire 
incremented to include the following five effects; (1) 

Hotspot profiles ace accounted for by using a pattern facto~ 
of 0.3 for the first stage stator of the high pressure 

turbine and a value of 0.13 for all other rows of coole^ 

airfoils in the turbine (s). (Pattern factor is defined as 

the ratio of the difference between the hotspot and the 
average row inlet temperature to the difference between the 


average row inlet temperature and the combustor inlet 
temperature.) (2) The correction to the hot gas temperature 
due to the dilution of upstream cooling air is obtained froa 
a mass averaged enthalpy from which a revised gas 
temperature is calculated. (3) Currently, the correction 
due to relative velocity is merely an approximation; the 
relative total temperature is assumed to equal 0.92 of the 
absolute total temperature. This value is a modification of 
the 0.90 value found in reference 3 which included dilution 
from the stator cooling air. (4) Because work is extracted 
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froD the ']ds streaa, the downstreai roiis of cooled vanes aad 
blades are subjected to a lower gas tenperat ure . (5) The 
gas temperature used to determine the required cooling flow 
is first increased by 150® B to provide a safety factor for 
the cooled turbine airfoils. 

Allowable bulk metal temperature depends upon the year 
of material technology, the desired life and the duty cycle. 
(Duty cycle includes mission, atmosphereic conditions and 
power settings.) In an attempt to simplify these 
relationships the following equations relate 
state-of-the-art bulk metal temperatures to both life and 
year of material technology. (See APPENDIX A for the 
definition of the variables.) 

TNVAMB = (10. *YEAB -17640.) - ( 1 00*log (ELITE) - 400. ) 

TNBLAD = (10. *yEAB3- 17740.) - ( 100*log ( ELI FE) - 400. ) 

These equations have been obtained from infornation 
similar to that contained in figure 8 of reference 5. In 
that figure metal tempecaturo is presented as a function of 
the year of material technology. The figure shows both a 
turbine blade material temperature of 2110® B for 1985 
technology and also a historic increase of almost 10® B per 
year in the the allowable metal temperature. Stator blades 
are herein assumed to operate at a temperature 100® B higher 
than the rotor blades so that they will operate at 2210® a 
for 1985 technology. 

The cooling air temperature reflects the stage of the 
compressor from which the bleed flow is eKtracted. 

Dimensi o nl e ss cooling f lo w. - After the relative 
cooling flow, FACTOH, has been chosen from table I for the 
row of airfoils, and after the cooling ef feet ive ness, PHI, 
has been t?valuated from the above temperatures, then the 
dimensionless cooling flow can be determined from the 
following equation. For example, for a FACTOH OP 1.0 and a 
cooling ef fectivenass of 0.5, the dimensionless cooling flow 
is 0.022 for a full coverage film cooled turbine blade. 

«(I)/WG = FACTOK*0.022* (PHI/(1.-Pai)) **1.25 

The form of this equation is derived from a heat 
balance across the surface of a turbine blade in which 
turbulent flow is assumed for both the hot gas and the 
cooling air. The constant 0.C22 is obtained by 
extrapolating experimental data presented in figure 23 of 
reference 6 for a full coverage film cooled turbine vane. 

To account for endwall, shroud and disk cooling, and 
Leakage, the dimensionless cooling flow for each row of 
airfoils is increased by i factor of 4/3. Although this 
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value will vary froa engine to engine, its use seems 
reasonable when comparisons are made with flowrates 
predicted by other methods. 

Cooled t urbine e fficiency. - First, the 

uncooled-stage efficiency, EFF2, is calculated from the 
unccoled-turbine efficiency, BPFl, using equation (6«5.6) of 
reference 7. This calculation assumes that all 
uncooied-stage efficiencies and all stage pressure ratios 
are equal. 

EFF2= ( 1 - (1-EFF1* (1-BP<'*GAn) )♦* (1/XN) )/(1-RP**(GAH/XM) ) 

next, the dimensionless cooling flows from the stator 
and the rotor, excluding the endwall, shroud, disK, and 
leakage flows, are then used to determine the decrease in 
the thermodynamic efficiency of the turbine stage. This 
decrease in stage efficiency is equal to the product of the 
dimensionless cooling flow, one of the last two columns in 
table I, and the uncooied-stage efficiency. Subtracting the 
change in stage efficiency from the uncooled-s* >ge 
efficiency, EFF2, will yield the cooled-stage efficit :y, 
EFF3. 

EPF3 = EFF2- (W (I) /WG) *DELV*BFF2- (K (1+1) /WG) *DEL!I*EPF2 


Finally, the thernodynamic cooled- turbine efficiency, 
EFF4, is reconstructed from the cooled-stage efficiencies 
using equation (6.5.5) also of reference 7. This equation 
allows a different value of cooled-stage efficiency for each 
stage. 

XN 

EFF4= (1-TT(1-EFF3 (I)* ( 1 -RP»* ( 1/XN ) ) **GA H) ) ) / ( 1-RP**GAM) 

1 = 1 

Flow Chart of the Algorithnj 

The sequence for calculating these quantities can be 
followed with the aid of APPENDIX A and APPENDIX B. The 
former is a symbol list and the latter is a listing of the 
new subroutine. 

As shown in APPE'IDIX 8, before the «D0 20" loop, the 
uncooled thermodynamic stage efficiency, EPF2, is calculated 
from the uncooled thermodynamic turbine efficiency, BPF1. 

The "DO 20" loop separates the kind of cooling for the 
turbine, KINDOF, into the cooling configuration for each row 
of the turbine, NPACT(I). For example, in a one stage 
turbine with a value of KINDOF egual to 86, the values of 
MFACT(I) equals 8 and NFACT(2) equals 6. This indicates 
that the stator is cooled by a full coverage film and the 
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rotor is cooled partially by fila and partially by advanced 
convection. See table I for NFACT(1)=8 and NPACT(2)=6. 



The "DO SO** loop deteraines for each row of airfoils 
the effective hot gas temperature, T, the cooling 
effectiveness, PrtI , the required cooling flow for each row, 
W(I), the change in stage efficiency, DELM, the accuaulated 
required cooling flow made dimensionless by the flow into 
the turbine, PCBLED, and the cooled-stage efficiency, EPP3. 

The "DO 90" loop reconstructs the cooling 
configuration, KINDOP, from NPACT(I) to reflect any changes 
made to NFACT(I) inside of subroutine COOLIT. Such changes 
occur when either a blade row requires cooling but was not 
originally intended to be cooled or when a blade row does 
not require cooling but was originally intended to be 
cooled. For the former case the blade row is assumed to 
have convection cooled airfoils; for the latter case the 
blade row is assumed to ha^e solid, uncooled airfoils. See 
table I where NPACT (I) equals 1 and 0 respectively. 

Finally, beyond the "DO 90" loop, the cooled-turbine 
efficiency, EFFh, is reconstructed from the cooled-stage 
ef f iciencies. 

Input to the Algorithm 

I nput to g ene r al t he rm od ynamic cycle c ode . - Input to 
the algorithm ace discussed below. The names marked with an 
asterisk are not always required because of the reasons 
cited. 

♦EPF1 is the theraod yna mic efficiency of the uncooled 
turbine. It is only required if the thermodynamic 
efficiency of the cooled turbine is to be calculated- If 
EFF1 is not known then it can be calculated by using the 
cooling conf igurations of the airfoils, the corresponding 
dimensionless cooling flows, and the above three equations 
relating EFFI, EFP2, EFF3, and EFF4. Or it can be found by 
a trial and error procedure by varying the input value of 
EFFI until the known output value of EFF4 is achieved. An 
example of such a calculation might be to determine the 
required cooling flow and cooled- turbine efficiency for a 
transpiration cooled turbine Knowing the corresponding 
values for a film cooled turbine. 

•ELIFE is the desired life of the turbine airfoil. It 
is assumed to be 10 000 hours if the design values of TflVANE 
and TMBLAD are known. 

♦FARCX is the fuel to air ratio of the cooling flow. 
It is used to evaluate the specific enthalpy, HC. If a gas 
property subroutine is not available and the diluted gas 
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nanufacturers as long as the proper cooling configuration 
and bulk netal temperatures are used. 

To determine the significance of the new cooling 
algorithm, a comparison was made between its predictions and 
those of the old version of the COOLIT subroutine for an 
advanced engine. The comparison of predicted cooling flows 
shews that the new cooling algorithm predicts a flow which 
is about 65 percent higher than that predicted by the old 
COOLIT subroutine- in?, major reason for the difference in 
the predicted flows is the fact that the old COOLIT 
subroutine did not allow for either a safety margin or a 
hotspot temperature. 

E ffect of coo l ing f low on p erformance . - In the new 
algorithm, compressor bleed flow can affect engine 
performance in two ways. (1) As the cooling flow exits fro* 
the turbine blades and enters the gas stream it causes a 
loss in total pressure of the mainstream due to increased 
drag and thereby causes a decrease in turbine efficiency. 
(2) Any cooling flow which is ejected into the gas stream 
downstream of the rotor is directly chargeable to the 
engine cycle since work cannot be extracted from it. 

The new cooling algorithm predicts a higher required 
cooling flow than did the old COOLIT subroutine- It 
therefore follows from the previous paragraph that NNEP in 
conjunction with the new cooling algorithm predicts lower 
thrust levels than that predicted by NNEP in conjunction 
with the old COOLIT subroutine. In one example the net 
thrust decreases by over 9. 2 percent because of the extra 
cooling flow which bypassed the rotor. 

Mhen NEWEFF-.TRUE. , the turbine efficiency is degraded 
by cooling f^cj injection. When NEMEFF-.PALSE. , turbine 
efficiency is not degraded by cooling flow injection. i';.e 
difference between these two options for the case cited 
above is a decrease in net thrust of 1.1 percent. Host of 
this difference is due to the decrease in the efficiency of 
the high pressure turbine (HPT) because most of the total 
cooling flow is required for the HPT. For the case cited 
the efficiency of the HPT decreases 3.76 points while the 
efficiency of the LPT decreases 0-88 points- 

It is interesting to differentiate between the relative 
error in thrust attributed to the decrease in turbine 
efficiency and that attributed to the cooling flow which 
bypasses the rotor. The error introduced into the thrust 
level due to inaccurate cooling flow is several tines the 
error due to uncorrected turbine efficiency. 
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APPENDIX A. - NABES OP VARIABLES IN SUBROUTINE COOLIT 

AIF 1 Air content of the hot gas at the entrance to a 

ro« of airfoils# Ibn/sec 

AXR2 dir content of the bleed flow froa the row of 

airfoils iaaediately upstreai, Iba/sec 
CALBLD indicator which is set to .TRUE, when bleed flow 
requirement is to be determined 

DSLN percent change in thermodynamic efficiencyof a cooled 

stage for each one percent of required cooling flow; 
change in thermodynamic stage efficiency due to the rotor 
DELV change in thermodynamic stage efficiency due to the stator 

DH power extracted from each turbine stage, Btu/sec 

EFPI thermodynamic efficiency of an uncooled turbine 

EFF2 thermodynamic efficiency of an uncooled stage of a turbine 

SFF3 thermodynamic efficiency of a cooled stage of a turbine 

EFP4 thermolynamic efficiency of a cooled turbine 

ELIFE desired life of the turbine airfoil 
FACTOR cooling flow weighting factor 
FAHCI fuel to air ratio of the cooling flow 

FARG fuel to air ratio of the hot gas stream at the 

inlet to a row of airfoils 

FARGX fuel to air ratio of the hot gas stream at the 
inlet to the turbine 

FUEL1 fuel content of the hot gas stream at the entrance 
to a row of airfoils, Ibm/sec 
FUEL2 fuel content of the bleed flow from the row of 
airfoils immediately upstream, Ibm/sec 
GAB (GAMAVG-1) /GABAVG 

GABAVG (GABIN+GAB0UT)/2 

CABIN specific heat ratio of the hot gas as function of 
TIN and FARGX 

GABOUT spev.ific heat ratio of the hot gas as function of 
TOUT and FARGX 

HC enthalpy of the cooling flow, Btu/lbm 

HG enthalpy of the hot gas stream at the inlet to a 

row of airfoils, Btu/lbm 

KGUP enthalpy of the hot gas stream at the inlet to a 

turbine, Btu/lbm 

HPT power extracted by the turbine, hp 

II number of rov s of cooled airfoils in a turbine 

which exceed nine 

KVANE indicator which determines whether a row of 
airfoils is v. stator or a rotor 
KINDOF an ordered combination of digits representing the 
cooling configurations throughout the 

the turbine starting with the stator of the first stage 
B BOWS 

NEHEPP indicator which is set equal to .TRUE, when a new 
thermodynamic efficiency for the cooled turbine 
is to be calculated 

NFACT(I) digit representing a particular cooling 

configuration for a particular row of airfoils in 




d cooled turbine 

N whole number of stages in the turbine 

N1STT component number of the first turbine 
PHI general name for PHIC(X) 

?HIC(X) cooling ef feet iveness evaluated at the local value 
of the effective gas stream temperature 
PC3LED accumulated cooling flow required to cool the 
turbine, made 

dimensionless by the flow entering the turbine 
PR pressure ratio across turbine, inlet/exit 

PROFIL combustor pattern factor; or radial 
temperature profile 

RP pressure ratio across turbine, exit/inlet 

ROWS number of rows of cooled airfoils in the turbine 

SAFETY safety margin of 150®R used when determining the 
required cooling flow 
STAGES number of stages in the turbine 
sntlBLD sum of dimensionless bleeds from all turbines; 

approaches unity as engine design converges 
T argument of function PHIC(X); 

average effective gas temperature,® B 
TAOT temperature ratio across turbine as calculated 

by stage efficiencies 

TC total temperature of the cooling flow at the point at 

which it is extracted from the compressor,® R 
TG gas temperature at the inlet to a row of airfoils,® R 

TIN total temperature of the gas into the turbine,® R 

TH general name for either TflVANE or TM3LAD,®B 

T.1BLAD allowable bulk metal temperature foe rotor blades,® R 

TBVANE allowable bulk, metal temperature for stator vanes, ®R 

TOUT total temperature of the gas out of the turbine,® R 

WCOOL dimensionless cooling flowrate for a row of cooled airfoils 

W(I) the cooling flowrate for a row airfoils, Ibm/sec; 

or the cooling flowrate 

for a row of airfoils in the turbine plus 
endwal 1, shroud , and disk cooling and leakage, Ibra/sec 
WG gas flowrate at the inlet to a row of airfoils, Ibm/sec 

WG1 WTFLO, Ibm/sec 

WTFLO total gas flow into the turbine, Ibm/sec 

XDELN general name for DELN 

XFACT general name for FACTOR 

XN N 

YEAH first year of service for stator vane material; 

called YEARV in INPUT SUBROUTINE 
YEARS first year of service for rotor blade material 



I ..1,1 




appendix B. 


listing CF SUBBOUTINS rA.Ou: 
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SOBRODTINE COOLIT /STAGPq ifTMn^ = 

'?^^°.f«HCX,r,aGX,PB| ' ““°’^'”*'T0»T.TC,PCBLED,EPPl,ErP4,HP 

implicit REAL*8 (A-Ii.o-Z) 

COMMON /JBLBBD/ SOMBLD YEAB Pi r»p 

COMMON /BLADE/ YEARS 

common /PP/ PB0FIL,N1STT * ^BLAD, NEHEFF 
logical CALBLD 
logical newepf 

m?Pfa2t"iS/o;5,L'o°’, s'fJ^'JOl .KOns ,20, 

= (*-Tn, /DABs'd- re,- 5/ 

N=AMAX0 (1,N) 

GAMIN=THEHM(5,TIN,FARGX) ®^OG10(ELIFE)-400.) 

GAMA7G- (GAMIN + GAMOOT) /2 

gam= (gamavg-d/gamavg 

TA0T=1 
RP=1/PS 
XN = N 

EFF2= (l.-(l. -EFF1*n -Rp**r»uA . 

DH=HPT/FL0AT(N)/1.4t; <^-/^N))/(^-'RP*MGAH/XN)) 

HC-THERU (4,TC,FARCX) 

SAFETY'150. 

FARG=FARGX 

tg=tin 

HG=THERM(4,TG,FARG) 

PCBLED=0- 
HG1=WTFL0 
WG = WG1 

hgop=hg 

kvane=-i 

R0»S = 2*STAGES ♦. 1 

m=rohs 

M=AMAX0(2,M) 

II=M-9 
DO 20 1=1, H 

00 TO 10 

kindof=ki^!);%”"""" (KINDOF/10) 

GO TO 20 
MPACT(M+1-I) = 1 
»(I) =0. 

DO 00 1=1, H 
KVAME=KVANB* (-1) 

COOLING FLOW WEIGHTING 'FACTOR* A<>«:Tr-*Pn 

rnv,ioK ASSIGNED TO *KINnOF' 


13 




i 



c 


c 

c 


30 

c 

c 

c 

c 

40 

c 

c 


c 

c 

c 

c 


5 0 

60 

65 


c 

c 


70 


COOLING CONFIGURATION 
N=NFACT(I) 

F\CTOR=XFACT(N) 

IF (KVANS,LT.O) GO TO 30 
TM=TMVANE 

DIHENSIONLESS CHANGE IN STAGS EFFICIENCY FOR A ONE 
PERCENT CHANGE IN STATOR COOLING FLOW 
DELN=XDSLN (N) 

GO TO 40 
TM=T?JBLAD 

DIMENSIONLESS CHANGE IN STAGS EFFICIENCY FOR A ONE 
PERCENT CHANGE IN ROTOR COOLING FLOW 
DELN=XDELN (N* 10) 

APPROXIMATION FOR RELATIVE TOTAL TEMPERATURE 
TG=-92*TG 

CORRECTION TO *T* FOR SAFETY MARGIN OF 150 DBG 

T=TG+SAFETY 

PHI=PHIC(T) 

CORRECTION TO *PHI» FOR PATTERN FACTOR/RADIAL PROFILE 
•PEOFIL* HAS VALUE OP 0.30 POE THE FIRST STATOR OF THE HPT 
PH 1= (PROFILE PHI) / (PROPIL+1.) 

PHI=DMAX1 (PHI,0.D0) 

IF THE BLADE IS MEANT TO BE COOLED AND IT NEEDS TO BE COOLED 

IF (NFACT (I) .GT.O.AND. PHI.GT.O.) GO TO 50 

IF THE BLADE DOES NOT NEED TO BE COOLED 

IF (PHI.EQ.O. ) GO TO 60 

IF THE BLADE NEEDS TO BE COOLED 

BUT WAS NOT MEANT TO BE COOLED 

FACTO R=2. 0 

NFACT (I) =1 

W (I)=WCOOL (PHI) *FACTOH*WG 
GO TO 65 
NFACT (I) =0 

DELN=DELN*W (I) /WG*EFF2 
IF {. NOT. NEWEFF) DELN=0 
IF (KVANE.GT.O) DELV=DELN 
CORRECTION FOR ENDWALL COOLING/LEAKAGE 
W (I) =W (I) *4/3 

NSW FOEL/AIR RATIO AFTER DILUTION FROM UPSTREAM ROW 
AIRl = WG/( 1.4-FARG) 

FUEL1=WG*FARG 
AIR2=W (I) / (1. +PARCX) 

FUEL2 = W (I) *FAHCX 

FARG= (FUEL 1* PUEL2) / (A I R 1 ♦ AIR 2) 

IF (KVANE.GT.O) GO TO 70 
EFP3=EPP2-DELN-DELV 

TA0T= (1.-EFF3* (1. - (BP**(1-/XN) ) ♦♦GAM) ) ♦TAUT 
HG=HG-DH/WG 

HG= (WG*HG*H(I) ♦HC) / (»G*W(I) ) 

TG=THEHM ( I^HG, FARG) 

WG = WG>W (I) 

PROPIL=, 1 3 

PCBLED=PCBLED*W (I) /WG 1 


1 4 




1 


I 






80 


90 












appendix C, - CHANGES TO SUBROUTINE TURBIN 

Aftt'r the iit ate sent 
•'INPLICIT... •• 

on line 2 , insert the stateaeot 
"DIflENSION KIND(20)**, 

2. After the statement 
"C0?*..10N /J3LEED/..." 

on line 10, insert the two statements 

^Con.-ION /BLADE/YEARV,rnVANE,TnBLAD,NEWEFF- 
"ConMON /PF/ PRUFIL,N ISTT". 

3. After the statement 
"LOGICAL CALBLO" 

on line 11, insert the statement 
"LOGICAL NEWEFF", 

After the statement 
"EFFT=DATINP (6, JCX) " 

on line 45, insert the two statements 

and 13ATOUT (8,JCX) =0ATINP(11,JCX) " 

"IF (NEWEFF.AND,NIT.EQ.0) DAT0UT(11,JCX)=DATINP(11,JCX)". 

Replace the statement 
"EFF=DATINP( 11,JCX) •' 
on line 95 with the statement 
"EFF=DATINP(8,JCX) 

6 . After the statement 
"TOTE.NF (JP1) =TOUT” 

on line 109, insert the six statements 
"IF (N13TT.S0-0) NISTT^JCX", 

"IF (JCX, E3. N ISTT) PROFIL=0.3« 

'Ir (HPT.EO.O) K IN D ( JCX) =D AT I NP { 1 4, JCX) " 

"KIND0F = KIND(JCX1", V . A) , 

"IF (HPT.EO.O) GO TO 12V* 
and 

"IF (KI N DOF, EO. 0.) KINDOF=83*». 

7. Delete the two statements 
"PACT0R=DATIN?( 14, JCX) '• 

on line 110 and 

"IP (FACTOR. EQ.O) FACTOR =1.0" 
on line 111. 

9, Change the term "NSTAGE" on line 113 to "STAGE". 

9. Change the terms "NSTAGE, FACTOR" on line 114 to "STAGE, KINDOF" 

10. Change the term "LED)" on line 115 to the terms 
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"LED, DATINP(1 1,JCX) , DATOOT ( 8 , JCX) , H PT, «T 1 , FA R 2 ,FA R 1 , PR) 

11. After the statenent 
"IP (DABS... JO TO 60" 

on line 121, insert the statenent 

"IF (DABS(BFP-DATOOT (8, JCX) ) . (IT. 1.P-4) JO TO 60", 

12. After the statement 
"SUf1BLD = S0HnLD..." 

on line 122, insert the statement 
"CONTINOB", 

13. After the statement 
"DATOUT (7,JCX) ..." 

on line 133, insert the statement 
"IF (IDONB(JCX) .EQ.O) GOTO 125", 

14. Change the statement 
"DATOtJT (9,JCX) =PR" 

on line 135, to 

"125 DATOUT (9, JCX)=PR". 

15. After the previous statement on line 135, insert the statement 
"DATISP (14, JCX) =KINDOF". 


I 
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APPENDIX D 


CHANGES TO SUBROUTINE INPRT 



r 

f 




! 

I 




1. After the statement 
"COnflOM /JBLSED/... '• 
on line 11, insert the statement 
"COMJIOM /BLADE/YEARS, TNVANE, THBLAD, 


NEWEPF'*, 


2. After the statement 
”10) I,DATINP(2, 1) " 
on line 98 insert the statement 
•'KINDOP = DATINP ( 14,1) 


3- After the term "DATIN? (2 , 1) ” on line 99, insert the 
terms •*,KINDOP, THVANE, TMBLAD”. 


4. Replace on the rigb t-hand-side parenthesis, ") ", with a hyphen 
and add to format statement 520 the continuation 
"1.* THE COOLING CONFIGUBATION IS ,I9,» TBTAMB IS 

•,F6.0,* TBBLAD IS SF'.O)”. 

(Note the blank after the parenthesis.) 


I 


1 a 



APPENDIX E.- CHANGES TO SDBROOTINB INPUT 
K After the term "CALBLD" OD line 16, insect the tern ",NBWEFF'*. 

2. After the statement 
"COHNON /JBLSEDA.." 

on line 23, insert the two statements 

"COnnON /BLADE/ YBARB, TH7ANE, TMBLAD, NEWEFF" 

and 

"COMnON /PF/ PROPIL ,N 1STT**. 

3. Replace the term ’♦YBAR'* on line 28 with the terms 
'•NEWEFP, YEARS, YBARV«. 

4. After the statement 
•’CALBID=.FALSB. •• 

on line 38, insert the three statements 
»»NBWEFF=.FALSE.*», 

••YEARV=1985.", 

and 

••YEAR B= 1985. •». 

5- Delete the statement 

"YEAR=1985." 

on line 39. 

After the statement 
"READ (8,D)« 
on line 50, 
insert the statement 
"N1STT=0’'- 

7. After t ne statement 
"READ (8, D) " 

on line 76, insert the statement 
"YEAR=YEAR7". 
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